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We report on an analytical description of spin-dependent electronic transition rates which are 
controlled by a radiation induced spin-Rabi oscillation of weakly spin-exchange and spin-dipolar 
coupled paramagnetic states (S = |). The oscillation components (the Fourier content) of the net 
transition rates within spin-pair ensembles are derived for randomly distributed spin resonances 
with account of a possible correlation between the two distributions that correspond to the two 
individual pair partners. The results presented here show that when electrically or optically detected 
Rabi spectroscopy is conducted under an increasing driving field B ± , the Rabi spectrum evolves from 
a single resonance peak at s = £l R , where Q R = 7B1 is the Rabi frequency (7 is the gyromagnetic 
ratio), to three peaks at s = £l R , s — 2Q R , and at low s <C Q R . The crossover between the two 
regimes takes place when Cl R exceeds the expectation value 5 of the difference of the Zeeman 
energies within the pairs, which corresponds to the broadening of the magnetic resonance lines in 
the presence of disorder caused by hyperfine field or distributions of Lande g-i actors. We capture 
this crossover by analytically calculating the shapes of all three peaks at arbitrary relation between 
Q R and S . When the peaks are well-developed their widths are as ~ 6%/Q R . 

PACS numbers: 42.50.Md,76.30.-v,71.35.Gg 



I. INTRODUCTION 



Over the past decade, pulsed electrically detected mag- 
netic resonance (pEDMR) spectroscopy has been used 
increasingly for the investigation of the physical nature 
of spin-dependent electronic transitions which influence 
conductivity such as excess charge carrier recombination 
or transport transitions through locali zed paramagnetic 
states as seen in amorphous inorganic^l, crystalline^^ 
as well as organic^^ semiconductor materials. In most 
of these experimental studies, pEDMR exper iments are 
conducted within a pulse-probe scheme^^ where the 
current through the host materials of the spin systems of 
interest is measured while short and intensive magnetic 
resonant pulses are imposed. The pulsed magnetic res- 
onant radiation prepares coherent spin non-eigenstatcs 
from initial, well defined eigenstates before the pulse. As 
the changed spin-states will also cause changes of the 
spin-dependent conductivity, one can gain information 
about the prepared coherent spin-state by integration of 
the electric current transient after the pulse which even- 
tually, on long time scales (compared to the length of the 
coherent excitation), will return to its pre-pulse steady 
state. The probed charge (obtained from the integrated 
current transient) depends on the coherent spin state af- 
ter the pulse which in turn depends on the pulse pa- 
rameters (length, frequency, intensity!^. Therefore, a 
measurement of the charge as a function of the applied 
pulse length can reveal the propagation of a spin state 
in presence of the resonant radiation pulse. Thus, cur- 




FIG. 1: Schematic illustration of the Rabi oscillations with 
frequency Q R in spin-| pair. Components of the pair, a and b 
have different environments causing random shifts, S a and S b , 
from the resonant frequency, ui . Relevant for PEDMR are 
the initial and final spin configurations | tt) an d | only. 

rent detectable observations of spin Rabi oscillations arc 
possible. 

Most EDMR detectable spin-dependent electronic 
transitions reported in the literature are due to Pauli- 
blockade effects which occur for transitions between two 
paramagnetic states with S = 5. These systems, illus- 
trated in Fig.[l] usually require weak spin-orbit coupling 
as found in materials with low atomic order numbers 
(this means silicon and carbon materials) as well as suffi- 
ciently weak spin-spin coupling (which means exchange 
and dipolar interaction) within the formed pairs. In order 
to allow one of the electrons within this pair of paramag- 
netic states to undergo a transition into the other para- 
magnetic state (thereby forming a singlet spin manifold 
due to the Pauli exclusion principle), the spin pair state 
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before the transition requires non-negligible singlet 
content (('015') ^ 0) for the transition to have a non- 
negligible probability. The special nature of such inter- 
mediate pair controlled spin-dependent transition rates 
was first recognized by Kaplan, Solomon and Motifs who 
explained the magnitude of continuous wave EDMR ex- 
periments at the time. With the advent of pEDMR about 
a decade ago, this model also became most significant for 
the understanding of many of the EDMR detected coher- 
ent spin motion experiments. 

When pEDMR is applied to the intermediate pair pro- 
cesses described by Kaplan, Solomon and Mott, the ob- 
servable applied to the spin ensemble is permutation sym- 
metry (the singlet operator), it is not the magnetic polar- 
ization of the spin ensemble as it is the case for conven- 
tional magnetic resonance spectroscopies which are based 
on the detection of radiation. Some implications of this 
observable change have been discussed theoretically in 
previous studies o n inte rmediate pairs for cases of no in- 
trapair interaction d 12 * 15 !, cases of weak but non-negligible 
exchange interaction^ (weak here means that the ex- 
change interaction is much smaller than the spin-Zeeman 
splitting of the pair partners but not necessarily weaker 
than the difference of the pair partners Larmor frequen- 
cies), and for cases where disorder within the spin ensem- 
ble^ is significant. These numerical studies have shown 
that an electrically detectable spin-Rabi oscillation can 
contain various harmonic components which can essen- 
tially form a "fingerprint" for the spin-Hamiltonian of 
the observed pairs. Thus, conducting a Fourier analysis 
of an observed spin-Rabi signal (one could call this Rabi 
spectroscopy) can give microscopic information about the 
nature of charge carrier states or of paramagnetic defects. 

Most of the previously published pEDMR stud- 
ies h ave been c onducted as Rabi spectroscopy experi- 
mentd 2 ^ 7 * 8 * 10 * 11 !. For most of these experimental data, 
a correct interpretation would be impossible without the 
in forma tion provided by the existing theoretical stud- 
i e ji2 | i5| tiT] Nevertheless, these studies can only provide 
limited support for an experimental analysis since nu- 
merical simulations can only provide answers about the 
behavior of a simulated system for a fixed set of param- 
eters, but they do not reveal analytical expressions that 
can be fit or directly compared to experimental data, 
and most importantly, they oftentimes do not enhance 
the qualitative understanding of a simulated system, ft 
has been the goal of this study to overcome this problem 
by finding a closed analytical form for the description of 
spin-Rabi oscillation controlled spin-dependent transi- 
tion rates within spin pairs with S = h. The expressions 
derived in the following reveal the dependence of all har- 
monic components found with electrically, and similarly, 
with intermediate spin-pair controlled optically detected 



transition rates on the parameters of both, the observed 
physical system as well as the performed Rabi-oscillation 
experiment. 

The parameters of a performed Rabi-oscillation ex- 
periment are the spin S = J Rabi oscillation frequency 
fl R = "fB\ which depends on the driving field strength 
Bi as well as the gyromagnetic ratio 7. The parameters 
characterizing the observed spin pair are given by the 
spin-orbit controlled respective (/-factors for each pair 
partner as well as local, and hyperfine fields which in 
general are different for the two pair partners as well. 
For our purposes, all these parameters can be taken into 
account by the difference S of the pair partners Larmor 
frequencies. Taking this into account, we antic ipate from 
the previously reported numerical simulation ;] 12 ! 15 ! that 
in the weak-driving regime where Q R ^ d, the Fourier 
transformation F(s) of a Rabi oscillation transient ex- 
hibits only one peak at s — Q, R , while in the strong- 
driving regime fl R 3> S) there is one peak at s = 2Q R 
and no peak at s = fl R . The crossover between these two 
situations occurs around Q R ~ i5. 

This crossover with increasing Q R was demonstrated 
experimentally for spin-dependent polaron pair recombi- 
nation pro cesses in different organic semiconductor ma- 
terials^^. Currently, it is rather difficult (if not impos- 
sible) to extract quantitative information from such ex- 
perimentally measured Rabi spectra F(s) due to the lack 
of theoretical predictions. On the other hand, the the- 
oretical problem is well-posed. When the exchange and 
dipolar coupling strength within a pair are smaller than 
5 (which we assume) , the shape of the Fourier transform 
depends only on two parameters: S , which is the r.m.s. 
value of 5, and Q R . 

In the following, the Fourier components of the Rabi 
spectrum F(s) are calculated analytically. We find that 
the width of the crossover region is broad and extends 
from ^ 0.3, where the s = 2fl R peak appears, to 

w 2, where it dominates over s = Q R peak. Most 
importantly, the analytical treatment reveals that F(s) 
consists not only of two peaks (as discussed in most ex- 
perimental studies) but rather of three peaks. The origin 
of the third peak, which occurs at frequencies s <C £l R , 
is due a disorder induced distribution of 5, which implies 
that even at strong f2 H , the two spins in the pair do not 
precess entirely in phase. This third peak is harder to 
observe experimentally compared to the two with higher 
frequencies, yet a prediction of its shape and evolution 
with £l R is provided. Also, a correlation study of disor- 
der within pairs affects the shape of an ensemble average 
of F(s). This follows from previous numerical study of 
the ratio between the magnitudes of the s = Q R and 
s = 2Q R Rabi oscillation peaks^ which was conducted 
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for particular value sets of S and £l R . 

The study presented here is outlined as follows. In 
Sect. II, a qualitative derivation of conductivity changes 
in pEDMR experiments is given which reproduces the 
result of rigorous consideration in Ref. In Sect. Ill, 
the analytical expression for a disorder averaged Fourier 
transform of A<t(t) is presented before this is used in 
Sect. IV for the consideration of correlation effects be- 
tween the two random distribution of the two pair part- 
ner resonances. An analysis and discussion of these re- 
sults is then presented in Sect. V. 



II. THE DEPENDENCE OF PEDMR INDUCED 
CONDUCTIVITY CHANGES ON THE PULSE 
DURATION t 

The pair partners within a spin pair are denoted with 
a and b respectively. Before a pulse is applied to a spin 
pair, it will rest in one of its four spin-cigcnstatcs, as 
both a and b can be either in a | \.) or in a | f) state. 
The important qualitative observation made in Ref. [T5] 
is that only initial configurations | 44) and | ft) 01 the P a i r 
exist as the other two eigenstates with singlet content are 
very short lived. The | 44) and | ft) states are therefore 
responsible for the change of conductivity after the end 
of the pulse. Without confinement of generality, one can 
assume that at t = both a and b are in the 1 4-) state. 
The respective populations of the | 4-) states will then 
evolve according to the Rabi formula 



n a ,b(t) 



-fa 
4V 



(1) 



where S a = uJ a — uj and 5b = ujf, — uj are, the detuning 
frequencies of a and 6, respectively. A detuning frequency 
is the difference of a Larmor frequency uj a or ujb and the 
frequency uj of the driving field. After the pulse ends at 
t = t the pair is in the | 44) state with probability P^ = 



n a (r)ni,(r) and in |ft) state with probability Pt 



n 



l- 



n a ("?")) (l — n b(T~))- Then the probability to find the pair 
in one of the states | 44) or | ft) is equal to 



P(r) = P n + P 



44 



1 - n„(r) - n b (r) + 2n a (r)n b (r). 



(2) 



It is easy to see that Eq. ^ also applies when the pair 
is initially in | ft) state. Eq. ^ coincides with the cor- 
responding expression for the r-dependent part of the 
diagonal elements of the density matrix of the pair de- 
rived in Ref.Q^l The probability P(t) serves as the initial 
condition for the transient restoration of the steady state 
current after the pulsed Thus, Act can be identified 



with P(t) within a factor. An expression for A<t(t) av- 
eraged over the contributions from all pairs within a pair 
ensemble is then obtained from 



(Aa(r)) 



1 



dS a d5 b exp 



2(5„ 2 



Aa(S a ,S b ,r). 

(3) 



Note that this expression for (A<t(t)) can also be written 
as 



(Aa(r)) = l-2P(r) + 2T 2 (r), 
where the function T(t) is defined as 

T(r) =^L- ( dS e- s2 ^ s o ( 



(4) 



3in 2 ^(n2+* 2 ) 
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In the limit of long pulses T{t) approaches a constant 
in an oscillatory fashion; the amplitude of the oscilla- 
tions falls off slowly, as r -1 / 2 , with the length of the 
pulsed For strong disorder (5 ^> f2 fl ) the derivative 
can be expressed through the zero-order Bessel functiorP 
T'(t) = 2- 3 /V/^lS~ 1 J (Q R r). 

While the second term in Eq. Q describes Rabi os- 
cillations within either component a or b of the pair, the 
third term "knows" about the collective spin precession 
of a and b. However, the T 2 -term also contains contri- 
butions from the individual precessions of pair partners 
a and b. We will therefore subtract these contributions 
and group them with the T-term in Eq. Q prior to per- 
forming the Fourier transform. By substituting Eq. ([T]) 
into Eq. |2]), we get 



Acr(r) = 



f 



G 1 (<y.,<y»,T) 



2 2(fia+£)(na +5 2) 
+ G 1 (6 b ,6 a ,T)] + G_(5 a ,5 b , t) + G + (S a ,S b , r), 



(6) 



where the functions describing the three harmonic Rabi 
oscillation peaks are defined as 



G 1 (S a ,S b ,T) = 



, (7) 



G_(<5 a A,r) 



fin 4 



cos 
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G + (5 a ,5 h , T ) = 



tin 4 



n 



% + S* + VS%+S*)t} 



(9) 

The above terms G 1: G + , G_ describe the peaks s = fl R , 
s = 2f2 H , and s <C fl R , contained in the Fourier transform 
F(s), respectively. 



III. AVERAGING OVER DISORDER WITHIN A 
SPIN PAIR ENSEMBLE 

Variations of the magnetic resonance frequency of each 
pair partner in each individual pair can occur due to: 
(i) variations of the spin-orbit coupling which changes 
the g-factoM This is seen in disordered materials 
where the lengths and angles of chemical bonds can vary 
strongly^ 2 -!, (ii) Due to random hyperfine fields which can 
strongly fluctuate throughout a material because of the 
small polarization nucl ear sp ins even at low temperature 
and high magnetic field d 10 * 2 ^. We define the Fourier spec- 
trum of the conductivity change from the steady state as 



F00 



<2tcos(st)((Act(t) - Acr(O)) 



The expression can be decomposed into three contribu- 
tions F(s) = Fj(s) + F (s) + F 2 (s), which derive from the 
terms G u G_, and G + terms in Eq. (|6j). Obviously, the 
time integration of each term yields a (5-function. 

Our task is to perform the averaging of each (5-function 
over disorder, as in Eq. pj). We start from Fi(s), which 
describes a peak near s = fl R . For this contribution, the 
averaging over 5 a , 5b reduces to the product of averages 



M*) 



8% 



5{s - ^fl + P) 



x / dS b e -^ 5 « 
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5 2 
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It is convenient to evaluate the integral over <5„ with the 
help of the (5-function. The second integral can be re- 
duced to the error-function leading to 
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exp 
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(12) 



where the function / is defined as 



- Vne b erfc(v^) 



(13) 



For Fa(s), it is convenient, after substituting ^ into 
( |10[ ), to perform an integration over S a , 5b in polar coor- 
dinates. Upon introducing the new variables 



5 a = \/2vcos(j), 5 b = v2wsin</) 
the expression for F2(s) acquires the form 
fit 



(14) 
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dvve' v /s ° 



d({> 



5[s - ytfll + 2v 2 COS 2 ( 



Q 2 + 2v 2 sin 2 



(Q.% + 2v 2 cos 2 <f>) (n% + 2v 2 sin 2 



(15) 



Without the denominator, it is straightforward to per- 
form an integration over (/>, which yields 

J d<p 5 (s - ^n 2 + 2v 2 cos 2 $-yJ n 2 + 2v 2 sin 2 <j^j 

(16) 



(^-v 2 f -W R s 2 {W R + V 2 ) 



(10) With the denominator, Eq. (16 1 is to be divided by the 



value \{s 2 — 2S1 2 — 2t> 2 ) 2 of the denominator where the 
argument of the (5-function is zero, which yields 
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5 2 



dvve y2 l s a 
\s 2 -2n 2 R -2v 2 \ 
1 



-v 2 ) -ills 2 



(17) 



Eq. (17) is defined only for s > 2i} R . As s approaches 



2Q R , both brackets under the square root turn to zero. 
At the same time, the integration interval also shrinks 



to zero. However, the factor Is 



20 2 



2v 2 



the 



denominator is nonsingular near s = 2Q R . To illuminate 
the behavior of F2 (s) near the threshold, it is convenient 



to make the substitution w = 



- Q 2 R in the integral 



of Eq. (17). It then assumes the form 



(f-n H ) 2 




dw e~ w/s o 



(f-fi*) 1 



\f-2w\J (% + n R y-w 



(18) 
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Now we see that only the first two factors in the denom- 
inator are singular when s is close to 2f2 R , where the 
s = 2Q R peak occurs. In this domain we can set w = 
in the last two factors and take them out of the inte- 
grand. Then, the remaining integral readily reduces to 
the modified Bessel function, I (y), and we get 



F,(«) = 



U 2 * 2 







[ e «8 ' 




(s - 2ft fl ) 2 
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(19) 



where 



G(6) 



da; e 



7re 



-6/2 r 

-'o 



(20) 



The analysis of the shape of the s = 2fl R peak given in 
Sect. V. will reveal that the approximation in Eq. (19) 
describes not only the vicinity (s — 2f2 R ) <C £l R but the 
entire peak when Q R is bigger than O.3<5 . 

Finally we turn our attention to the peak at s <C £l R - 
The initial expression for F 2 (s) differs from Eq. (15) 
for F (s) only in one respect. Instead of the sum, 
■J £1 2 + 5 2 + + S 2 , in the argument of the (5-function 
it contains a difference, \J Q R + 5% — \/Q R + §f. The an- 
gular integration in polar coordinates is therefore per- 
formed in a similar way as Eq. ( 16 1 



><)[s - + 2v 2 COS 2 . 

\2Sll + 2v 2 - 



Q 2 + 2v 2 sin 2 . 



,2V 



fi 2 



(ni + v 2 ) 



(21) 



Note that, compared to Eq. (21), the integral Eq. (16) 
has an extra factor of 2. This is because the argument 
of ^-function in Eq. 
(j), while in Eq 



(21) turns to zero at two values of 
(Il6| , it turns to zero at four values of (f>. 
To get the final expression for F„ one has again to divide 
Eq. (21 ) by the value of denominator at the zeros points 
of thed-function, which is equal to \(s 2 — 2fl 2 R — 2v 2 ) 2 , 
i.e., the same as in F 2 (s). This leads to 



Fo(s) 



fl% / dvve~ v2 / 6 o 

^ Jy/($-*) 2 -ni* 

V J T+sn R 



\s 2 -2n 2 R -2v 2 \J{Q 2 +v 2 ) 



(22) 



We see again that only the first factor in the denominator 
is singular at the lower limit v — (%■ + sf^) 1 ^ 2 . Thus, 
for small s, the peak is described by substituting s = 
into the second and third factors in denominator, which 
leads to 



Or 

4<J„ 2 



2 /r2 

dvve~ v '°° 



« J^r R ^(f-v 2 f 



(23) 
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The remaining integral can be expressed via the Macdon- 
ald function, K (y), which yields 



7Mg K 



SQ R 

S 2 

"(J 



(24) 



We will see in Sect. V. that the expression in Eq. (24) 



describes the entire peak when fl R is big enough, fl R > 6 



IV. CORRELATION BETWEEN THE PAIR 
PARTNER DISORDER 

Due to the proximity of the pair partners in each pair, 
the disorder related randomness of S a and S b may be cor- 
related. Examples for such a correlation could be the 
common exposure of polaronic state in organic semicon- 
ductors to an overlapping nuclear spin bath or the corre- 
lation of the spin-orbit interaction in a disordered semi- 
conductor due to local strain fieldfP. In such cases, o„ 
and S b are not statistically independent and the degree 
of overlap can be expressed by a correlation parameter, 
x (0 < x < 1). Then the joint distribution function of S a , 
S h assumes the form 



2ir5 2 Vl 



: exp 



2x5 2 8\ 



25 2 {l~x 2 ) 



(25) 

We study the effect of correlation for the limit Sl R S 
when all three peaks are well-developed and do not over- 
lap. In this limit, we can use the expansion 



n n + in the arguments of the (5-functions. We can 
also replace w^ 2 + S 2 b by Q R in the denominators of 

Eqs. ([7]) to ([9]). With these simplifications the expres- 
sion for F x (s) assumes the form 



Fi(s) 



1 



^ini^/i^x 2 

S 2 + 6 2 

x exp 



2S 2 (1- x 2 



d5 b j d5J 2 b 

2x6 „S „ . 




2fT R 



(26) 



G 



Integrating over 5 a with the help of the ^-function yields 
F,(.s) = ; / (!>'' 



exp 


Sl R (s-S~l R ) 






(40AVi - x 2 } 


^/2fi R (s-Q R ) 



x 5 2 exp 



s 2 b 

252(1 - x 2) 



cosh 



5 b x^/2Q, R (s - Q, R ) 
S 2 (l~x 2 ) 



(27) 



Subsequent integration over S b is straightforward leading 
to 



Ms) 



4Q Ry /n R (s - n R ) 



2x 2 n R (s - n R ) 



x exp 



ft 



We see that the correlation parameter, x, enters this ex- 
pression only in the prefactor. In the limit of x 0, Eq. 



(28) matches Eq. (12) as can be seen when Q R ^> S is 
assumed and Eq. (12) is expanded around s = fl R . 
In the limit of large fl R the definition of F 2 (s) becomes 

F 2 (s) = ( — 1 ! / ,M ,M 



165 Vl - a: 2 

S 2 + 8 2 - 2xS b S, 



x exp 



2fi H 



(29) 



Both integrals over 5 a and <5 t can be taken explicitly upon 
introduction of polar coordinates 



r = \/2(5 2 + S 2 ), <p = arctan f 



s a -s b 

S a + 6 b 



(30) 



so that Eq. ( 29 1 assumes the form 
1 



Ms) 



32<S Vl - .t 

2 



<ir r < 



x exp 



r 2 (l + x cos 20) 
4<5 2 (l-a; 2 ) 



8 [s - 2n R 



r 



(31) 



After integrating over r by using the <5-function, the re- 
maining integral over (j> reduces to I (y), and we arrive 
at 



Ms) 



ttCI r 



exp 



n R (s - 2n R ) 
' <5 2 (i-^ 2 ) 



X In 



Q R x(s — 2Q R ) 
S 2 (l-x 2 ) 



(32) 



Similarly to the non-correlated case, F 2 (s) is expressed 
through I a {y). However, note that the argument of I {y) 
in Eq. (|32| is completely different from Eq. (fl9|. In fact, 



Eq. d32h was derived for the case when I (I) inEq. (|20j 
should be replaced by 1. 

In the presence of disorder correlation, the shape of 
F (s) can be also expressed via the Macdonald function 
with x-dependent argument. To see this, we take the 
definition 



Ms) = 



x exp 



1 

166 Vl - x 2 
S 2 + S 2 - 2x8 b 5, 
25 2 {\-x 2 ) 



dS„ dS b 



5 s- 



2n R 



(33) 



(28) and introduce polar coordinates 



S 2 + 8 2 - 2x6 b S a 
V 15 2 {l-x 2 ) 

arctan 



l + x 



1 -x 



(34) 



Eq. ( 33 ) then becomes 
1 



Ms) 



drre 



d<p S I s ■ 



n R 



sin 2(f) 



The integration over <f> can be done explicitly by using 
the 5-function, yielding 



Ms) 




dr 1 



1/2 / 6*(1- X 2) 



(36) 



The reduction to the Macdonald function can then be 



achieved by the substitution r 2 
veals 



Ms) 



n R 



ss 2 Vi 



K 



S 2 VT 



x* 



r x which re- 



(37) 



V. ANALYSIS AND DISCUSSION 

The most important results of this study are the ana- 
lytical expressions for the lineshapes of Rabi-oscillation 
peaks as they can be found in the Fourier analysis (the 
Rabi spectra) of the pulse length dependent conductivity 
changes Act(t) that are measured with pEDMR experi- 
ments. For the case of uncorrelated disorder, Eqs. (12), 
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(191, and (24), reveal these peak shapes for the oscilla- 
tion peaks Fjfs), F (s), and F 3 (s), respectively. For the 
case of correlated disorder, the same peaks are described 
by Eqs. (28), (32), and (37). We consider now the limit 
Q R ^> S , when all three peaks are well developed. It is 
easy to see that, for uncorrelated disorder, all three peaks 
exhibit the same exponential tail 



Ms) 



Fo(s) 



7T 1 / 2 

F 2 (s) 



7T 



8\/2sas 8as 
where the characteristic width of the tail is given by 



(38) 



AS 



51 



(39) 



Naturally, all three peaks shrink with increasing il R . It is 
less trivial to realize that, for the same deviation from the 
origin, the peaks F , and F 2 have a larger magnitude than 
Pj, whose terms contain Q R 3> as in the denominator. 
This relation between the peaks is illustrated in Fig. [2];. 

We consider now the opposite limit of strong but un- 
correlated disorder where 8 3> fi s and x = 0. Eq. 
(19) implies that F 2 (2f2 R ) = is always finite at its 

peak. At the same time, the value of F 1 (s) at s = 2fl R 
is , i.e., it is bigger than F 2 (2Q R ). Therefore, the 

s = 2Q R peak is indistinguishable on the background of 
the s = Sl R peak. This behavior reflects the physics of 
the weak driving regime^ 2 -' where only one component of 
the pair can be in resonance with the driving field at 
any time. Numerically, however, the s = 20 R peak is 
pronounced already at Q R > 0.3# , as shown in Fig. [2k. 
For the given strength of Sl R , the approximation Eq. ( jl9[ ) 
using the modified Bessel function is already justified. 

The low frequency Rabi-oscillation peak which is de- 
scribed by Eq. p4| diverges logarithmically in the limit of 
s — >■ when F (s) oc ln(l/s). Nevertheless, the peak still 
loses to Fi(s) in the weak-driving regime due to its small 
prefactor . The lineshape of this peak is described by 



Eq. (24) when the ratio Q R /8o exceeds 1. This could be 
the reason why this peak has not received much atten- 
tion in previous experimental studies since this regime is 
hard (yet not impossible) to attain experimentally. 

Fig. [3]illustrates how the intrapair correlation of disor- 
der affects the shapes of the peaks in the strong-driving 
regime. One can see that the prime effect of correla- 
tion is a dramatic narrowing of the low-frequency peak. 
This narrowing reflects the fact that the low-frequency 
peak is entirely due to inequivalences of the pair partners, 
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FIG. 2: The shapes of the three peaks s <g £l R (red), s = 
fl R (blue), and s — 2Q. R (green) in the Fourier transform of 
Acr(r) are plotted for three values of the dimensionless Rabi 
frequency w = Q. R /8 : (a) w = 0.3, (b) w = 1.1, (c) w = 2. 



which are suppressed by the correlation. Another effect 
of strong correlation is that F x (s) develops a maximum at 
(s — Q R ) = as/2 and a minimum at (s — f2 R ) = (1 — x)as. 
The origin of the maximum is that, for full correlation 
(x = 1) the portion of resonant pairs in which only 
one partner participates to the Rabi oscillations van- 



ishes. More precisely, F 1 (s) oc (s — fi Ji ) 1 / 2 exp 
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FIG. 3: Illustration of the effect of intrapair correlation of 
disorder on the shapes of the three peaks s <C £l R (red), s = 
Q R (blue), and s — 2£l R (green), in the Fourier transform of 
Act(t). The dashed lines correspond to uncorrelated disorder 
(x = 0), while the full lines correspond to strongly correlated 
disorder (x = 0.95). The dimensionless Rabi frequency is 
w = tt R /S = 1.2. 



at x = 1, giving rise to a maximum. For small but finite 
(1 — x) such pairs exist and cause "normal" divergent 
behavior F 1 (s) oc }~ x at s — >■ f2 B , resulting in a min- 



imum. Finally, the s — 2Q R peak becomes enhanced by 
correlation near the origin for the same reason why F^s) 
gets depleted. Conservation of the total area, which ap- 
plies for all three peaks, is achieved due to depletion in 
the body. 

We note that disorder broadens the peaks in F(s) only 
to the right (the higher frequencies) from corresponding 
thresholds. This is due to the adopted definition [Eq. 
(10)] of the Fourier transform 

definition, [F 2 (s) + F 2 (s)] 1/2 , 



. If we used the standard 
where F(s) is defined as 



F( S ) 



drsin(sr)((Acr(r) - Acr(0)) 



(40) 



the disorder will broaden the peaks both to the left and 
to the right (to higher and lower frequencies) from the 
threshold. This is illustrated in Fig. [4] for the central 
peak at s = Q R . A formal expression 
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2tt 51 s 



dS 



S 2 - ft 2 - 



5 2 J \25% 



(41) 



emerges as an obvious modification of Eq. (12 1. This 



integral in Eq. ( 41 ) can be expressed through the error- 
function. It is nonzero for s < Q R and for s > £l R . Near 
,s = fl R it diverges as \s— f2 n | -1 / 2 . As shown in Fig.JZJ the 
inclusion of F(s) into the Fourier transform (as done for 
most previously published pEDMR and pODMR Rabi 
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FIG. 4: Solid lines: the Fourier transform of the central 
peak plotted using the definition [F 2 (s) + F 2 (s)] 1/2 from Eqs. 
(10 1, (401 for two values of the dimensionless Rabi frequency 
w = Q. R /5 . The dashed lines show the corresponding Fourier 
transform plotted using F(s) only. 



spectroscopy studies) leads to a significant unnecessary 
increase of the peak width which complicates the analysis 
of experimental data. 

For the illustration of the artificial broadening ef- 
fect, we refer to an experimental data set which is dis- 
played in the inset of Fig. [5| The data shows electri- 
cally detected spin-Rabi oscillation in a 7r-conjugated 
polymer diode measured at room temperatures on iden- 
tically prepared samples and conditions as for the ex- 
periments described in Ref. 23] The diode consisted of 
an ITO/PEDOT/MEH-PPV/Ca/Al device stack. The 
signal was measured by detection of changes to a for- 
ward steady state current of / = 100^A current, under 
application of a 4V bias at room temperature. The ex- 
periment was conducted under an applied magnetic field 
of 344mT, the applied magnetic resonant excitation had 
a frequency of 9. 6606GHz and a driving field strength 
Bi = 0.519(36)mT produced by a 27W coherent pulsed 
microwave source. The inset of Fig. [5] shows the mea- 
sured raw data. It displays a rapidly dephasing oscilla- 
tion which is due to inhomogencitics of the applied B\ 
field, Ref. The main plot of Fig.[5]shows two different 
Fourier transformations of the data displayed in the inset 
plotted on a normalized scale as a function of frequency 
in units of fl R = r yB 1 . The black plot displays the abso- 
lute Fourier transform of the experimental data while the 
red plot displays one the in-phase (cos) component of the 
Fourier transform. From these plots, it becomes evident 
that the absolute Fourier transform displays significant 
broadening of the Rabi spectrum without providing ad- 
ditional insight. It shall be noted that the peak widths of 
the plot in Fig. [5] are too close to the frequency resolution 
of the Fourier transform to unambiguously identify that 
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FIG. 5: Fourier transformation of an electrically detected 
spin-Rabi oscillation measured in an organic polymer diode. 
The inset displays the measured data set, for details about 
the measurement see text as well as Ref. [23] The black data 
of the main panel displays is the absolute Fourier transform 



F(s) 2 + F(s) 2 of the data in the inset. The red data dis- 
plays is the real part F(s) of the Fourier transform of the data 
in the inset. The blue line represents a Fit of the experimental 
data with the analytical function F(s) = Fo(s) + Fi(s) + -^(s) 



the peaks of the red plot are more asymmetric compared 
to the black plot. With the given length of the original 
data set (500ns), an additional symmetric broadening ef- 
fect is introduced to both Fourier transforms which over- 
laps the intrinsic peak shapes. However, for the data sets 
displayed in Fig. [5] one can conclude, that the analysis of 
experimental spin-Rabi spectra measured with pEDMR 
and pODMR should be conducted on Fourier data ob- 



tained from Eq. ( 10 ) 



In order to scrutinize the analytical description of the 
Rabi-spectrum of a pair process given above, we can now 
fit the experimentally obtained Fourier transform of the 
electrically detected electron spin-pair Rabi-oscillation 
displayed in Fig. [5] with the analytical expressions in 
Eqs. ( 12 ), ( 19 ), and (f24j) for the Fourier transforms of an 



uncorrelated spin pair ensemble Pi(s), ^(s), and Fq(s), 
respectively. It is known that weakly coupled polaron 
pairs in MEH-PPV do not exhibit a correlation of their 
local hyperfine fields^ so an agreement between the ex- 
perimental data set and the sum of Eqs. (12 1, (19), and 



(24) is expected. Since the expressions for Fj(s) are ab- 
solute, the fit of F(s) = F (s) + Fi(s) + F 2 (s) has only 
three fit parameters, namely (i) a global scaling factor 
of the unitless ordinate, (ii) the disorder factor f2 R /<5„ as 
well as a third fit variable which is the width of a ho- 
mogenous (Lorentzian) broadening of F(s). This broad- 
ening is due to an experimental artifact, coming from a 
Fourier transformation of a finite time interval as well 
as the measured decay of the Rabi-oscillation which is 
known to be due to inhomogeneities of the experimen- 
tally generated B\ fielcP^. Thus, the length (500ns) of 
the experimentally recorded Rabi oscillation determine 
the frequency resolution of the resulting Fourier trans- 
form while the decay (« 150ns) of the Rabi-oscillation 
further contributed to homogeneous broadening of the 
experimental Fourier transformed data. 

Using the three fit variables for F(s) mentioned 
above, we obtain a good fit of the experimental data 
as displayed by the blue line in Fig. [5j In this fit, 
the disorder parameter 5 /£l R = 1.54 with an esti- 
mated error of less than 0.2. Given the Rabi nutation 
frequency Vt B = r yB 1 = 14.5(1. 0)MHz, we obtain 
S = 22.3(3.3)MHz. This value is equivalent to a hyper- 
fine broadened polaron line width of $0/7 = 0.80(14)mT, 
a value that is in good agreement with the widths 
of EDMR d etected polaron spin resonances in MEH- 
pp^fHEEIM] We conclude from this that the analytical 



expressions Eqs. (12), (19), and (24) capture adequately 



the Rabi oscillations in weakly coupled spin pairs and 
should therefore be used for the fit of pEDMR data 
governed by non-correlated spin pair ensembles. 
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